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Abstract: We compute a one-loop effective action for the constant modes of the 
scalars and the Polyakov loop matrix of M = 4 SYM on S 3 at finite temperature and 
weak 't Hooft coupling. Above a critical temperature, the effective potential develops 
new unstable directions accompanied by new saddle points which only preserve an 
5*0(5) subgroup of the SO (6) global R-symmetry. We identify this phenomenon as 
the weak coupling version of the well known Gregory-Laflamme localization instabil- 
ity in the gravity dual of the strongly coupled field theory: The small AdSs black hole 
when viewed as a ten dimensional, asymptotically AdS5 x S 5 solution smeared on the 
S 5 is unstable to localization on S 5 . Our effective potential, in a specific Lorentzian 
continuation, can provide a qualitative holographic description of the decay of the 
"topological black hole" into the AdS bubble of nothing. 



1. Introduction 



The finite temperature behaviour of SU(N) gauge theories at large N provides a re- 
markable window into the physics of black holes and stringy gravity via the AdS / CFT 
correspondence Features of the phase structure of large N, H = 4 supersymmet- 
ric Yang-Mills (SYM) theory on S 3 at finite temperature are now known to mirror 
aspects of semiclassical gravity on asymptotically AdSs spacetimes |2|, [| [|, ||, §]. 
Importantly, the work of |J has demonstrated that studying the field theory in the 
tractable regime of weak 't Hooft coupling may allow us to deduce qualitative physics 
of the string theory dual emerging at intermediate and strong couplings. The weakly 
coupled field theory exists in one of two thermodynamically stable phases separated 
by a first order deconfinement transition; in addition there is the possibility of a 
thermodynamically unstable saddle point || [7|, § . Indeed, the strong coupling grav- 
ity dual also exhibits a first order Hawking-Page transition M between two stable 
geometries: thermal AdS space and the big AdS-Schwarzschild black hole, mediated 
by the thermodynamically unstable Euclidean small AdS black hole bounce. 1 



The aim of this paper is to take the qualitative matching of thermodynamic phase 
structure one step further and to explore certain questions which have dynamical 
consequences. Specifically, the small AdSs Schwarzschild black hole, when viewed 
as a ten dimensional asymptotically AdSs x S 5 solution smeared uniformly on the 
S 5 , has a classical dynamical instability (in addition to having a thermodynamic 
instability). This is a Gregory-Laflamme instability to localization [13|, [14|, [15], |16[| , 
wherein the small black hole originally smeared on the S 5 develops an instability as 
its horizon size in AdS 5 decreases below a critical radius (or equivalently, above a 
critical temperature). The unstable mode leads to the small black hole becoming 
non-uniform and eventually point-like on the S 5 , breaking the associated 50(6) 
isometry down to S0(5). The question we aim to answer in this paper is how, if at 
all, this dynamical instability to localization may be seen within the framework of the 
holographically dual thermal field theory on a three-sphere. Remarkably, within the 
regime of validity of perturbation theory we find that above a critical temperature 
at which the radius of the three-sphere becomes comparable to the Debye screening 
length, the weakly coupled field theory exhibits a clear signal of such an instability. 2 



The Af = 4 theory on S 3 at finite temperature has two scales, namely the radius 
R of the three-sphere and the temperature T. Consequently the theory possesses two 

1 The comparison between the phase structure of the weakly interacting thermal gauge theory 
and strong coupling gravity dual has recently been extended to include non-zero chemical potentials 
for the global U(l) 3 C S0(6) R charges in the field theory [|o) (see also Jl|]). The latter is dual to 
charged black hole geometries in AdSs. 



2 A similar phenomenon has been studied in 1+1 dimensional Yang-Mills theory in 17 
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tunable dimensionless parameters, the combination (TR) and the 't Hooft coupling 
A = gy M N. On S* 3 x S 1 , at weak 't Hooft coupling A C 1, we calculate a one 
loop quantum effective action as a function of homogeneous background expectation 
values for the Polyakov loop 3 and the six scalar fields transforming in the adjoint 
representation of SU(N). More precisely, we compute a finite temperature effective 
action of the Coleman- Weinberg type on a slice of the full configuration space pa- 
rameterized by the N eigenvalues of each of these fields. The one loop computation 
is valid for a wide range of temperatures at weak 't Hooft coupling 

< TR < - . (1.1) 
A 

A crucial point, which ensures the validity of perturbation theory at high temper- 
atures within the range above, is the infrared cutoff provided by the finite size 
of the three-sphere. In particular, we will mainly be interested in temperatures 
(TR) ~ 1/y/X which is well within the above range. At these temperatures the 
radius of S 3 becomes comparable to the Debye screening length (V\T)-\ In flat 
space, at these scales one needs to resum the one loop thermal mass of the modes to 
cure infrared divergences in perturbation theory. On S 3 however, the presence of an 
explicit infrared cutoff, namely R ~ (VA~T) _1 , ensures the validity of perturbation 
theory at these length scales, albeit in the parameter \/\ as opposed to A. Our cal- 



culation differs from the work of (10| in that we have generic values for the Polyakov 
loop matrix and scalar fields, but vanishing chemical potential. Thus we explore the 
landscape of the effective potential away from local minima. 

We now summarize the main features of the quantum effective action we obtain in 
equations ( |3.23|) , ( |3.24j ) and (|3.25|) of the paper. Firstly, the effective action provides 



a combined static effective potential on the space of the eigenvalues of the Polyakov 
loop matrix and of the constant modes of the scalar fields in the M = 4 multiplet. It is 
obtained by integrating out all fluctuations of these matrices (including off-diagonal 
ones) and all non-zero modes on S 3 x S 1 . For vanishing scalar fields, our action 
reduces to the unitary matrix model obtained in ||. Hence at low temperatures 
T <C R^ 1 , the effective potential exhibits a single saddle point with vanishing scalar 
fields and uniformly distributed Polyakov loop eigenvalues. This is the "thermal 
AdS" saddle point. At the Hagedorn temperature T# = — (R ln(7 — 4y / 3))~ 1 (or 
slightly below it for finite coupling) a first order deconfinement transition occurs 
beyond which a new global minimum appears with a gapped distribution for the 
Polyakov loop eigenvalues and vanishing scalar fields. This is the "big AdS black 
hole" minimum. 



3 We use the term Polyakov loop loosely since we are actually referring to the holonomy matrix 
of the gauge held around the thermal circle and not just its trace. 
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At temperatures above Th, the thermal AdS saddle point with vanishing scalar 
fields persists as a thermodynamically unstable saddle point, while the globally stable 
big black hole dominates the canonical ensemble. For all temperatures T# < T < 
R^ 1 /yf\, this continues to be the picture with the scalar fields forced to vanish by the 
tree level quadratic mass term arising from the conformal coupling to the curvature 
of S 3 . 

At temperatures T > R -1 / V\ however, something interesting begins to happen. 
The one loop contribution becomes comparable to the tree level mass while, crucially, 
the higher loop corrections remain perturbatively small, suppressed by powers of y/X. 
At a critical temperature T c ~ R^/y/X, a new unstable mode, along the R-charged 
scalar directions, emerges at the thermodynamically unstable "thermal AdS" saddle 
point. In the large N limit this critical temperature is found to be 

r « = ^ fl " 1 - (L2) 



Beyond this temperature, new unstable saddle points appear with lower action than 
the "thermal AdS" configuration which becomes unstable to rolling down to these 
new extrema. Importantly, at these extrema the scalar fields of the M = 4 theory 
have non-zero values. In the large N theory these expectation values are left invariant 
only by an 5*0(5) subgroup of the full SO (6) R-symmetry of the theory. However, 
the "big black hole" configuration continues to be the global minimum of the action 
with vanishing VEVs for all scalars which is consistent with expectations from the 
bulk gravitational physics. 




Small Black hole existence 



X ► 

Figure 1: Qualitative plot of small black hole (SBH) existence lines, and extrapolation of weak 
coupling critical temperature for the Gregory-Laflamme localization instability as a function of 
't Hooft coupling A. SBH comes into existence above the dashed line and below the Hagedorn 
temperature line. Localization occurs above the solid line. 
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We identify this high temperature phenomenon as a continuation to weak 't 
Hooft coupling of the Gregory-Laflamme instability encountered in the gravity dual 
of the strongly coupled gauge theory. There too, a thermodynamically unstable 
saddle point, namely the small AdS black hole develops a new dynamical instability 
to localization on the S 5 breaking the 5*0(6) isometry to SO (5). The main difference 
is that at weak coupling, at temperatures of 0(R~ 1 any "small black hole" type 
unstable saddle point disappears pfl or merges with the "thermal AdS" configuration 



In fact, as the temperature is increased, we find more and more unstable modes 
in the large N field theory around the thermal AdS saddle point. Whenever the 
temperature hits the critical values 

t? = (i.3) 



labeled by the positive integers /, a new unstable mode appears. This phenomenon is 
consistent with the observations of in the context of the instabilities of the small 



black hole where, as the temperature is increased, new tachyonic modes emerge. 
Importantly, the unstable modes involve only the homogeneous fields on S 3 , all 
Kaluza-Klein harmonics remain massive. 



A naive extrapolation of our weak coupling result suggests that the critical tem- 
perature for the onset of the instability decreases with increasing A. It is therefore 
plausible that for large couplings, the localization instability kicks in at a tempera- 
ture where the small black hole saddle point still exists and one obtains a Gregory- 
Laflamme instability. This is illustrated in Figure |l|. We must also bear in mind that 
weak 't Hooft coupling translates to string scale curvatures in the dual geometry and 
at these scales it is possible that the small black hole makes a Horowitz-Polchinski 
transition to a highly excited state of strings [T8|, ^ [Tj| before the onset of any 
localization instabilities. 



We should point out that the large N limit plays an important role in all of 
the above, since only in this limit can we legitimately speak of "SO (6) breaking 
expectation values" in the field theory formulated on the compact space. On a 
compact space, in the quantum theory we must integrate over all points in field 
space which are related by the global symmetry. 4 At large N the integration 
measure for this averaging procedure yields a contribution which is sub-leading in 
N. We discuss this in more detail later in Section 3.3.2. 

4 For a discussion of these issues particularly from the viewpoint of the small black hole gravity 
solution, see |L9j , 



- 4 - 



The outline of our paper is as follows. In Section 2 we review the general story of 
black hole instabilities in AdS space. Section 3 is devoted to a detailed calculation of 
the one loop effective potential in field theory and determining its saddle points as a 
function of temperature. We present our conclusions, interpretations and questions 
for future study in Section 4. The analysis of unstable directions of the effective 
potential is presented in an Appendix. 



2. Instabilities and AdS Schwarzschild Black Holes 



We begin by reviewing the physics of Schwarzschild black holes in asymptotically AdS 
spaces |9|, ||, ||. Of particular interest is the the five dimensional AdS-Schwarzschild 
black hole which is an asymptotically AdSs solution to the vacuum Einstein's equa- 
tion with a negative cosmological constant: 

ds 2 = -V(r)dt 2 + ^ + r 2 dnl (2.1) 

where V(r) = 1 + -^- — ^-(1 + |^-) and R is the radius of AdSs. The black hole horizon 
is at r = rh where V(r) vanishes. This black hole emits Hawking radiation. To find 
the associated Hawking temperature at which the black hole is in equilibrium with a 
thermal bath at that temperature, we employ the usual trick of going to the Euclidean 
section and requiring that there are no conical singularities. Taking t — > —ix, the 
resulting Euclidean metric is 

dr^ 

ds 2 = V{r)d X 2 + ^- ) +r 2 dnl (2.2) 

Requiring that the metric has no conical singularity at r = leads to the x coordi- 
nate being periodic with period (3 = fp^jp • Asymptotically, as r — ► oo, the space is 
Euclidean AdSs with the Euclidean time direction being a circle of circumference (3. 
This implies an equilibrium Hawking temperature 

l_ 2rl + R 2 
~ (3 ~ 2irr h R? ' {2J) 

Notice that for a certain range of T, there are two solutions for given by: 



txTR 2 . . _ . 

There are no solutions unless T > ^| = T , which implies that the AdS Schwarzschild 
black holes only exist for temperatures above this critical temperature. When T > T , 
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RT 




Figure 2: Temperature as a function of the horizon radius. 



there are two geometries corresponding to the two values of rh, the small and the 
large AdS black holes. The specific heat C v ~ -^p- — 2-kR 2 t'^ ^ negative for the small 
black hole and positive for the large black hole. This implies that the small black hole 
is in an unstable equilibrium with the radiation bath. At equilibrium, the Hawking 
radiation emitted and the radiation absorbed from the heat bath are equal. If the 
black hole emits a little more than it absorbs, it decreases in size, thereby increasing 
its temperature (because of the negative specific heat) which means it emits even 
more, implying the existence of a thermodynamic instability to decay to thermal 
AdS. If on the other hand, the small black hole emits infinitesimally less than it ab- 
sorbs, it will grow in size until it becomes the large black hole. The large black hole is 
in stable thermal equilibrium by virtue of its positive specific heat. In the Euclidean 
setup, the thermodynamic instability of the small black hole manifests itself in the 
existence of a non-conformal negative mode in the small fluctuation analysis. 



Although the small AdSs Schwarzschild black hole is thermodynamically unsta- 
ble, it suffers from no classical dynamical instability, in the sense that there are no 
tachyonic small fluctuation modes. This changes when we add extra dimensions along 
which the black hole is uniformly smeared. For example, in the context of type IIB 
string theory on asymptotically AdSs x S 5 background, the five dimensional small 
AdS black hole solution is smeared uniformly on the S 5 . In what follows, we will 
review how this 10 dimensional solution suffers from a classical dynamical instability 
to localization on the S 5 [16] . 



2.1 Gregory-Laflamme instabilities and the Gubser-Mitra conjecture 



When black holes are smeared along some extra non compact "internal" direction, 
often the onset of a thermodynamic instability is accompanied by a classical dynam- 
ical instability signaled by the existence of tachyonic modes in a Lorentzian small 
fluctuation analysis. In fact, this dynamical instability is of a Gregory-Laflamme type 



12]: The instability is to localization in the extra directions. We will now review this 



in more detail. Such a link between thermodynamical and dynamical instabilities 
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was conjectured by Gubser and Mitra [T3|, [L4[]. Further evidence for this connection 
was provided by Reall |15|] and we briefly review his argument by considering a black 



string in a five dimensional asymptotically flat space [TL6| 



The black string in five dimensions is translationally invariant in one spatial 
direction. It can be viewed as a four dimensional Schwarzschild black hole solution 
smeared uniformly along one extra non-compact direction. If we consider just the 
four dimensional black hole in the Euclidean section, the Euclidean Lichnerowicz 
operator has a negative eigen-mode 

Ag } V = W> (2.5) 
where i] is a real. The association of this Euclidean negative mode with the ther- 



modynamic instability signaled by a negative specific heat was made in If we 



consider static fluctuations in Lorentzian signature, the equations obtained are the 
same as the Euclidean ones: = A^ when acting on static fluctuations. 



However, ( |2.5| ) is not the relevant equation in the Lorentzian space. The lin- 
earized equation governing the small fluctuations about a Lorentzian 4D black hole 
solution is 

Ai 4) V = 0. (2.6) 

If, on the other hand, we are considering a five-dimensional black string and fluc- 
tuations which carry momentum along the translationally invariant fifth direction, 
njiv(x, z) = h ( jfu(x)e tkz , the relevant fluctuation equation involves a five-dimensional 
Lichnerowicz operator: 

^h^ = (A^+^)e^ = 0. (2.7) 

This is the same as ( p.5|) with k 2 = rj 2 . For this value of k, there is a "modulus" corre- 
sponding to this mode, which signals the onset of instability in Lorentzian signature. 



This was called the threshold unstable mode in [p"5| 



The above analysis may seem to imply that every small AdS black hole, when 
viewed as solution in AdSs x S 5 which is smeared on the S 5 should exhibit a Gregory- 
Laflamme like instability. However, the situation is more subtle if the internal 



smeared directions are compact, as was discussed in jig] . This is simply because 
the momentum in the compact internal dimensions is quantized and this leads to a 
constraint on the negative Euclidean eigenvalue for the existence of a threshold un- 
stable mode. For the smeared AdSs x S 5 Schwarzschild case the negative eigenvalue 
has to be ~ l< jp^ which are the eigenvalues of the Laplacian on the S 5 . This in turn 
imposes a constraint on the horizon radius of the Gregory-Laflamme unstable black 
holes. For example, only for = OAR does the I = 1 mode on the S 5 becomes 
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unstable [16]. This is smaller than r h — which is the largest possible horizon 
radius for the small black hole. In other words, although the small black hole exists 
for T > T = ^ = it only becomes Gregory-Laflamme unstable to localization 
on the S 5 for T > ^f-. 



2.2 Hawking-Page transition and the big AdS Schwarzschild black hole 



The AdS/CFT correspondence relates the bulk string theory partition function on 
asymptotically AdSs x S 5 geometries to the partition function of the boundary SYM 
theory. In the semiclassical supergravity limit, the bulk partition function gets con- 
tributions from saddle points which are classical solutions to the equations of motion. 
Depending on the boundary geometry, there can be more than one bulk saddle points 
and in such a situation a careful sum [g, [| over all the saddle points is required. 5 

For example, the appropriate boundary geometry to calculate the canonical en- 
semble partition function is S 3 x S 1 . For temperatures T > T = there are 
three bulk saddle points with this boundary behavior. These are the two AdS 
Schwarzschild black holes (the small and the large) and the thermal AdS geome- 
try, which is Euclidean AdS with a periodic time direction. For temperatures lower 
than the critical temperature To, thermal AdS is the only saddle point. The different 
saddle points correspond to different phases, and the locally stable saddle point with 
the least action dominates the canonical ensemble. 



The actions of the three saddle points are formally infinite because of the infi- 
nite volume of thermal AdS and AdS Schwarzschild black holes. However, it turns 




Figure 3: The actions for the big and the small black holes as a function of temperatures. 



out that the differences between the actions of the different solutions are finite and 
subtracting the thermal AdS action from the action of the AdS Schwarzschild black 
holes || yields finite results. Taking the thermal AdS action to be zero, the action 

5 For a review of AdS/CFT correspondence with an emphasis over the sum over different bulk 
saddle points, see BjJ ]. 
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for the AdS Schwarzschild black holes is given by 

,2 
h 
,2 



Vo\(S 3 )r 3 R 2 -r; 



SBR ~ '7(5/ h (^ + 2rV (2 - 8) 

It is easy to see that the action of the small black hole is always greater than both 
the action of the large black hole and that of thermal AdS (Figure |3|), so the small 
black hole never dominates the canonical ensemble. The action of the large black 
hole becomes less than that of thermal AdS at the Hawking-Page transition temper- 
ature T HP = At lower temperatures, thermal AdS dominates while at higher 
temperatures the large black hole dominates the canonical ensemble. 



3. Weakly- Coupled Af = 4 SYM at Finite Temperature 

Recently, weakly coupled SU(N) Yang-Mills theories on S 3 have been studied for 
large iV at finite temperature || [7j , where it was shown that these theories exhibit 
an interesting phase structure at weak coupling. The order parameter distinguishing 
these phases is the expectation value of the Polyakov loop, which is non-zero in the 
high temperature phase and is zero in the low temperature phase. 

In fact, the theories have an interesting phase structure even in the limit of 
zero coupling |J, ||. This is because if we define the free theory as a limit of the 
interacting theory, we need to impose the Gauss's law constraint which enforces the 
zero charge condition on a spatially compact manifold. On a finite S 3 , this effectively 
introduces interactions which are non-trivial enough to lead to a first order Hagedorn 
type deconfinement phase transition at high temperatures. This Hagedorn phase 
transition is believed to be the zero coupling analog of the Hawking-Page transition 
in the dual geometry at strong coupling. At non-zero but weak 't Hooft couplings it 
is expected that there will be two phase transitions, a deconfining transition followed 
by the Hagedorn phase transition. 

In the above phenomena at zero coupling were deduced from a Wilsonian 
effective action for the Polyakov loop matrix obtained by integrating out all non- 
zero modes on S 3 x S 1 . The effective action for the Polyakov loop matrix turns 
out to be a unitary matrix model. In the large- N limit it is natural to introduce a 
continuous spectral density function p(a) for eigenvalues of the Polyakov loop, which 
are distributed on the unit circle in the complex plane. It was shown in || that in the 
low temperature phase, the uniform distribution of eigenvalues of the Polyakov loop 
is a stationary point of the effective action and in fact is the minimum. The uniform 
distribution signals a confined phase of the field theory and is the weak coupling 
continuation of the thermal AdS dual geometry. 
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In the high temperature phase, although the uniform distribution remains a 
saddle point, it is not a minimum. The temperature at which the uniform distribution 
saddle point ceases to be a minimum is the Hagedorn temperature, which corresponds 
to the temperature of the confinement /deconfmement phase transition in the free field 
limit. Above the Hagedorn temperature, the minimum action eigenvalue distribution 
becomes nonuniform and in particular, is gapped, i.e. the spectral density p(a) 
necessarily vanishes on a subset of the circle. This implies a finite free energy for 
colored external sources and is naturally associated to a deconfined phase. This is 
the weak coupling version of the big AdS black hole. 

The small AdS black hole unstable saddle point cannot, however, be seen in the 
free field limit. It is expected to exist at weak, non-zero 't Hooft coupling || [7|, ||. 



3.1 Gregory-Laflamme in field theory 

Our goal in this section is to see the analog of the Gregory-Laflamme instability in 
the weakly coupled J\f = 4 SYM at finite temperature on an S 3 . As we discussed 
in the previous section, this instability is associated with the localization of the 
smeared small black hole on the S 5 , breaking the associated SO (6) isometry. In the 
dual field theory we would expect this phenomenon to manifest itself as an unstable 
saddle point which breaks the SO(6) R-symmetry. Such symmetry breaking has to 
be a subtle dynamical phenomenon for two reasons: i) The tree-level scalar potential 
on S 3 x S 1 precludes a non-zero classical VEV for the scalar fields in the M = 4 
multiplet; ii) field theories on compact spaces do not usually exhibit spontaneous 
symmetry breaking. 

To see the above phenomenon in weakly coupled M = 4 SYM we need to compute 
a quantum effective potential for the spatial zero modes of the scalar fields which 
are charged under the 5*0(6) R-symmetry. Thus in addition to the Polyakov loop, 
we will turn on the zero momentum modes of the six scalar fields and obtain a joint 
effective action for these degrees of freedom in perturbation theory. In addition, we 
will also see that spontaneous breaking of the global R-symmetry can occur on a 
compact space, due to the large N limit. 

Finally, it is important to note that we will compute the quantum effective 
potential along a special subspace in the space of field configurations, namely one 
where all the scalar fields and the Polyakov loop are simultaneously diagonalizable. 
This will be sufficient to explore the onset of the localization phenomenon which we 
are interested in. 
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3.1.1 The quantum effective potential 



In the canonical ensemble, the thermal partition function of a quantum field theory 
is equal to the Euclidean path integral of the theory with a periodic time direction 
of period f3 — ^ with anti-periodic boundary conditions for the fermions along the 
temporal circle. Since we are considering the Af = 4 SYM on S 3 at finite temperature, 
we perform a Euclidean path integral for the theory on S 3 x S 1 . 

We now describe the calculation of this path integral and hence the effective 
free energy (which is just log of the partition function) to one loop order. This will 
be a Euclidean lPI-effective action for static, spatially homogeneous fields and will 
therefore be interpreted as a static effective potential. We note here that although, 
strictly speaking, the static effective potential is a useful tool for identifying the 
equilibrium and locally stable field configurations, we will actually be interested in 
exploring its features away from local minima. 

Since the M = 4 theory is conformally invariant, the ratio of the radii of the S 3 
and the S l is the only physically relevant dimensionless parameter (in addition to 
the 't Hooft coupling constant) 

IH=* r - (31) 

where the radius of the thermal circle is f3, the inverse temperature and that of the 
S 3 is R. 

To do the effective potential calculation, it is necessary to expand all fields onto 
modes corresponding to the spherical harmonics on the S 3 and Fourier modes on the 
S 1 . We then compute an effective potential for the constant modes on the S* 3 x S* 1 by 
integrating out all other modes in the non-trivial background of the constant modes. 
In particular, we will compute the effective potential for the modes a and ip a , 

a = [ A), (3.2) 

and 

= , 3 . / 0a, a = l---6. (3.3) 
Vol(^) J 5 3 x5 i 

These are the spatially homogeneous, time independent pieces of A (the gauge field 
component along S* 1 ) and <f> a , respectively, and so are constant on S 3 x S 1 . 

Note that at the classical level, the Wilson line a is a true zero mode in the sense 
that the quadratic tree level action is independent of it. In contrast, the scalar modes 
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(p a are not zero modes of the classical theory since conformally coupled scalars have a 
mass of the order of the inverse radius of 5 3 , due to the curvature of the three-sphere. 

Tree level potential: 



The action for the bosonic fields is 

S h = [ d 4 x^Tr UF^F^ + \D^ a Dy a + ^<p a <Pa - \[<f>a, fa?) ■ (3.4) 

9ym J L J 

Here, summation over 50(6) indices is implied. The mass term above arises from the 

conformal coupling of the scalar fields to the curvature of 5 3 and in our conventions 

this mass is R~ 2 . 

In order to integrate out the non-zero momentum modes, we first shift the scalar 
fields by a background homogeneous mode 

4>a -> <Pa + 4>a ■ (3.5) 

The resulting theory has the following tree level effective action for the scalar modes 

5(0) = 0JivJ_ Tr (-{a^a? - h<p a ,<Ph] 2 + <?* Br 2 ) ■ (3-6) 
9ym \ 1 J 

In flat space in the absence of curvature induced mass terms, this scalar potential 

would have led to a Coulomb branch of vacua parametrized by the eigenvalues of 
a mutually commuting set of a and (p a . However, the conformal mass term for the 
{(fid} forces these to vanish at tree level since the potential (|3.6|) is a sum of positive 
definite quantities. In addition to the above terms there will be radiative corrections 
at finite temperature, obtained by integrating out the matter and gauge field fluctu- 
ations at one loop level. The nature of these radiative corrections and the question 
of stability of the ip a = classical solution will be the subject of the calculation below. 

Radiative correction at one loop: 



Computing the quantum corrections to the classical potential for generic homoge- 
neous background fields is technically difficult. Instead, we will compute the quantum 
corrected potential for mutually commuting background field configurations, namely 
those satisfying 

[a,<pa\ = [<fa,<fb\ = 0. (3.7) 
This will allow us to explore only a slice of the full configuration space. However it 
will be sufficient to address the issue of the presence of instabilities of the kind we 
are interested in. The classical potential on the space of these configurations is 

5<°) = [3Rti 2 ^— Tripl (3.8) 
9ym 
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Hence, as indicated above, for the effective potential calculation we shift 4> a — > 
f a + (pa where we assume that {<p a } and the Polyakov loop a are simultaneously 
diagonal. The fluctuations {4> a } which will be integrated out are fluctuations about 
the background, comprising of all the non-constant modes on S 3 x S 1 and off-diagonal 
components of the homogeneous modes. 

We denote indices along S 3 as i = 1,2,3 (not to be confused with the gauge 
index) and that along S 1 as 0. Unlike || who consider the case with vanishing 
background fields, we find it easier to work in a covariant gauge. In particular, we 
choose a conventional U R^ gauge" of a spontaneously broken gauge theory. To this 
end, we add to the action the gauge fixing term 

Sgf = [ dfix^gTrfoiA* + D A - i& a (j)a)\ (3.9) 

9ym J v J 

which implements a covariant gauge condition. Here and in what follows D = d +ia 
includes a, the zero mode part of Aq only. Also, we leave adjoint action by the 
background fields <p and a as implicit, i.e. <p4> = [ip, 0], <p 2 4> = [up, [ip, 0]], etc. 

We will now expand the gauge fixed action including the ghosts to quadratic order 
in fluctuations about the background fields (p a and the Polyakov loop a. The one 
loop determinants obtained by integrating out these fluctuations yield the radiative 
corrections to the effective potential. In order to proceed, it is particularly convenient 
to choose the Feynman gauge £ = 1. In this case the action for the bosonic quadratic 
fluctuations, including ghosts, takes a simple form: 



Sb 



9ym 



d 4 Xy /gTr \U (-D* - A« + <p 2 a )Ao + \H~Dl - A (l,) + 



\U-Dl - A« + R- 2 + rfWa + c(-D% - AM + 



(3.10) 



Here, A^^ and are the Laplacians on S 3 for scalar and gauge fields, respectively. 
The scalar Laplacian is simply VjV J acting on scalar functions whilst the vector 
Laplacian is defined via the quadratic terms in the action: 

J d^Xy/g Tr \F^ V F^ V = (3.11) 

J d 4 x^Ti( - lA^Dl + A^)# " U^oAo + V l A i ) 2 + 0(A 3 )). 

The second term on the right hand side is canceled by terms in the gauge fixing 
action (|3.9|) . From (|3.11|) , it follows that 

A (v) A i = Vj V j A 4 - tfjA* , (3.12) 
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where Rij is the Ricci tensor of S 3 . For the temporal component of the gauge field, 
the Ricci tensor does not contribute and the vector Laplacian is equivalent to the 
scalar Laplacian: 

A {v) A = A (s) A . (3.13) 

The eigenvectors of the scalar Laplacian are spherical harmonics Yg labeled by angular 
momentum £ G Z > with 

A^Y e = -£(£ + 2) R- 2 Y e , (3.14) 

and degeneracy (£+ l) 2 . The eigenvectors of the vector Laplacian can be split into 
two sets. Firstly, those in the image of V*, i.e. V l Yp with £ > 0, which satisfy 

A^V% = -£{£ + 2) R- 2 V% , (3.15) 

and secondly those in kernel, VjV^ = 0, also labeled by the angular momentum 
£ > 0, with 

A^V; = -{£ + l) 2 R- 2 V$, (3.16) 

and degeneracy 2£(£ + 2). 



We are now in position to integrate out all the bosonic fluctuations. For the 
vector modes, we must first split the A 1 into the image and the kernel of V*. To this 
end we write A 1 = B l + C l with Vj-B 1 = and C % = V*/. Integrating these out along 
with Aq, the ghosts and the six scalar fluctuations yields the following contribution 
to the one loop effective potential from bosonic radiative corrections: 

Sf = ilogdet £>0 (- J D 2 - AW + ifl) + ilogdet, >0 (- J D 2 - A« + cp 2 ) + 



+ (|-1) logdet,> (- J D 2 - + Vl) + 6 x ilogdet,> (- J D 2 - A (s) + R~ 2 + vl) ■ 

(3.17) 

We now describe the various contributions in some detail. The first two terms arise 
from integrating out the vector fluctuations Bi and Cj respectively where the sub- 
scripts remind us to exclude the £ = mode of A' 8 ' as discussed above equations 
fl3.15p and ( |3.16p . The third term in the effective potential results from the integrals 



over the ghost and Aq fields, each of which contribute the same factor with weight 
— 1 and | respectively. Finally we have the fluctuation determinants of the six scalar 
fields of the M = 4 theory. It is vital to realize that the determinants for the ghosts, 
the scalars and Aq all include the £ = modes as indicated. Note that the fluctua- 
tion determinants for the Cj, the ghosts and A combine neatly due to a complete 
cancellation between non-zero momentum modes yielding a simplified expression for 
the bosonic one-loop correction: 
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=\ logdet £>0 (- J D 2 - AW +<p 2 a )-± logdet £=0 (- J D 2 + y? 2 ) 

(3.18) 

+ 6 x ilogdet,> (- J D 2 - A« + R~ 2 + <p 2 a ) . 

We can make these formal expressions explicit by rewriting each of the terms above 
as a trace over the thermal frequencies or Matsubara modes and the discrete wave 
numbers on S* 3 . The result of the Matsubara trace of a typical contribution in the 
effective action takes the form (up to field independent additive constants) 

Trlog(-Z^ + e 2 ) = Tt(/3e-2^-e - ^cos(n/ta)) . (3.19) 

n=l 

Here e is the energy of the fluctuation in question with x° interpreted as imaginary 
time and we have used the fact that the operator Dq has eigenvalues (27rn/[3 + a)i; 
n G Z. The remaining trace on the right hand side is to be taken over the modes 
on S 3 and the gauge group. The gauge trace can be made explicit by using the fact 
that a and ip were chosen to be simultaneously diagonal 

a = /TMiag (6>i,0 2 , . . .6^), y? = diag (<p u <p 2 , . . . tp N ) , (3.20) 
which in turn yields in general 

N 

Tr f(<p, a) = f(<Pij> <Pij = ( fi- <Pj , <*ij = ^(Oi - 0j). (3.21) 

ij=l 

For SU (N) gauge group the 0i and ipi must satisfy 

N N 

^2(fi = 0, ^2 6i = mod 2ir. (3.22) 

i=i i=i 

These conditions are a consequence of the Hermiticity of the scalar and gauge fields. 
The 9i are each thermal Wilson lines for the Cartan components of the gauge field. 
Hence they can be shifted by an integer multiple of 2n by performing a topologically 
non-trivial gauge transformation which is single- valued in SU(N)/1*n. This must 
however be an invariance of the theory since there are no fields charged under the 
center Z N of SU(N). This results in the 6i being defined only up to an integer 
multiple of 27r. 

Including the explicit sums over the angular momenta on S 3 with the appropriate 
degeneracies, we can now write down the complete one-loop effective action resulting 
from integrating out all the bosonic fluctuations, 
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^ = E[4V^ + E^~ n ^ cos (™%) 

ij=l n=l 

+ m + 2) (f^+l) 2 i2- 2 +^- - E \e-N^^^U cos(^.)) 



n=l 

00 



+ 6 + I) 2 (f V /(^+l) 2 i?- 2 +^-- E ^e^v^W^J cos(n%) ) 

£=0 n=l 

(3.23) 



Finally, to complete the calculation of the effective action we have to include integrals 
over the four Weyl fermions. The effect of the background fields is very simple on 
these modes, it simply induces a mass squared (p 2 a for the fermions via their Yukawa 
couplings. The fermions are eigenf unctions of the spinor Laplacian on S* 3 which 
are also labeled by the angular momentum t > with eigenvalue — (£ + |) 2 and 
degeneracy 2£(£ + 1). The key difference between fermions and bosons at finite 
temperature is of course that the former obey anti-periodic boundary conditions 
around the thermal circle. Thus, when acting on fermions the Dq operator has 
eigenvalues (2-7r(n + 1/2)/ (5 + a)i; jigZ. 

Hence, with the anti-periodic boundary conditions around S 1 the fermionic con- 
tribution to the effective action is 



s i 1] = 4 E E 2 ^ + !) ( - f ^z+\) 2 R- 2 +vi 3 



ij=i e=i 

OO „ s n+1 



n=l 



(3.24) 



The complete one loop effective potential for the eigenvalues of the Polyakov 
loop variable a and those of the constant modes ip a is given by 



SribaiA] = S^+S^ + S^, (3.25) 

where the tree level potential for the eigenvalues is 

1 N 

S^=PR^^- E^L (3-26) 
9ym i=1 
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Table 1: The fields, their angular momenta, energy, degeneracy and their weight in the 
effective action. 



while and S^ 1 ' are determined in (3.23) and ( |3.24j ). Note the power of \j g\ M in 
front of the tree level term indicating that our one loop contributions are paramet- 
rically down by a power of the weak gauge coupling. 



Table 1 summarizes the properties of all the modes including their range of an- 
gular momenta, energies, degeneracies and the weight of the corresponding "log det" 
terms in the effective potential. 



Discussion 



We now make the following observations about our result: 



• Firstly, when the scalar fields <p a vanish we reproduce the results of 0. To see 
how this works, we note that the first term in the one-loop effective potential 
in the absence of the background scalars reduces to 

log II sin 2 (|) , (3.27) 

i<j 

which is precisely the Jacobian that converts the integration measure over 
the Hermitian matrix a into the appropriate measure for the Unitary matrix 

U = e ipa , 

J [da] det(-Do) = J dU = j J]^ J] sin 2 {^pj . (3.28) 

i=l i<j 

Here we have obtained this result following a different route from that adopted 
by |J. So when the scalars vanish it is easy to see that we reproduce the 
effective action written down in IjJ. 
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• Our effective potential was obtained by integrating out not only the non-zero 
momentum modes on S 3 x S 1 , but also the off-diagonal fluctuations of a and 
of the zero modes (n = t — 0) of <p a . (In principle one also integrates out 
fluctuations of the diagonal pieces but these give no nontrivial contributions, 
due to the adjoint action of the background fields on them). Note that the off- 
diagonal zero modes have masses \J R~ 2 + ^ a (<Aw — faj) 2 which are generically 
larger than the tree level masses for the diagonal components of the constant 
modes. Accounting for these constant modes is precisely what gives rise to 
the appropriate Jacobian factor for unitary matrices discussed above and the 
repulsive Vandermonde interaction between the eigenvalues (p a i of the scalars. 
The presence of the logarithmic repulsive force between the N eigenvalues (p a i 
is easily inferred from the fluctuation determinant for the i = n = modes in 

(El)- 

• The potential is only a function of the variable 

6 

¥&j = I>>»-^-) 2 - (3-29) 

0=1 

A non-zero expectation value for this will single out a particular direction in 
the space spanned by the six scalars and hence will spontaneously break the 
5*0(6) R-symmetry to SO (5) in the large N theory. Note that the notion of 
symmetry breaking on a compact space is unusual and we clarify this carefully 
in subsection 3.3.2. 

• The thermal interpretation of the various terms in the effective potential is 
also fairly clear. All the finite temperature contributions are accompanied by 
thermal Boltzmann suppression factors. In the zero temperature limit (5 — oo, 
these vanish leaving behind the terms which are manifestly linear in (3. The 
latter are the zero temperature or "Casimir energy" contributions on S 3 in the 
presence of scalar expectation values. We will discuss these zero temperature 
contributions in more detail in the following section. For now, we make one 
additional observation that the term proportional to \J<£ 2 a ij actually cancels 
against an identical piece emerging from the infinite sums over £. 

• The bosonic and fermionic determinants come with opposite signs as they must. 
It is, however, important to realize that they do not cancel against each other 
due both to finite temperature effects and the fact that the theory is formu- 
lated on a spatial three-sphere. A complete cancellation between bosonic and 
fermionic fluctuations only occurs in flat space and at zero temperature. We 
will come back to this point in the following subsection. 
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• The effective potential that we have obtained can be interpreted as a in- 
effective action for the eigenvalues of a and the <p a and yields the one loop 
partition function about a given configuration of eigenvalues. At local minima 
of the action, we may interpret the partition function as a thermodynamic 
partition function whose logarithm is the thermodynamic free energy. For 
generic points in the configuration space of field eigenvalues, the system will 
not be in a stable or static configuration. Nevertheless, we may formally define 
a free energy as 



We are interested both in extrema of this effective potential and the existence 
of and emergence of any new unstable directions as a function of temperature. 

• Finally, saddle points computed using the effective potential will be true saddle 
points in the full configurations space, even though we are setting to zero the 
background values of the off-diagonal constant modes. The reason for this is 
that, in our chosen background of diagonal fields, the off-diagonal fluctuations 
only appear quadratically and there are no terms linear in these fluctuations. 
This means that there are extrema of the full effective action where these can 
be consistently set to zero and these are the saddle points that we will find. 
However, there may well be other extrema where the fields p and a are not 
simultaneously diagonal. 



3.2 The zero temperature limit 

In this section we discuss the detailed form of the one loop effective potential at zero 
temperature. This is not relevant for the high temperature analysis which is our main 
focus. However, the limit of zero temperature will provide checks of our calculation 
and also important intuition about the behaviour of the one-loop corrections. 

The terms that survive in the zero temperature limit, j3 — oo, are basically 
Casimir contributions in the presence of background expectation values, and are 
identified as the terms without any Boltzmann suppression factors in the effective 
potential. The total Casimir free energy is therefore 



F(v,9) = S eS /{3. 



(3.30) 



N 




(3.31) 



ij=l 
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where 



CV) = ^^ + 2) v /(£ + l) 2 J R- 2 + ^ + 3 J2(£+l)W^+^) 2 R~ 2 + V 2 



1=1 e=o 

oo 



(3.32) 

Note the complete absence of any dependence on the Polyakov loop variable a. 

Furthermore, as we remarked earlier, the bosonic and fermionic determinants 
do not cancel against each other even at zero temperature though the system is 
super symmetric. This is purely a consequence of formulating the theory on a curved 
spatial manifold, namely an S 3 . Hence there is a nontrivial quantum correction to 
scalar potential of the M = 4 theory on S* 3 x R even though all supersymmetries are 
unbroken. The cancellation is recovered in the flat space limit. 

Evidently the infinite sums, at least individually, require regularization. This 
can be done by introducing a generalized ^-function, or simple Epstein series, 

oo 

S{z,9; V R) = [{n + 0/2n) 2 + { V R) 2 y z . (3.33) 

n=— oo 

In terms of this, the Casimir energy can be neatly expressed as 



C{^) = 2R~ L £(-§,0;pi2) -£(-§, 7r;pi2) 



~ + i) (£(-l 0; <pR) - S(-h vr; ipR)) 

(3.34) 

The bosonic and fermionic contributions mirror each other with opposing signs, the 
former yielding Epstein series with 9 = while the latter have 9 = it. 

It is easy to extract the zero point energy at tp — formally in terms of the 
Riemann zeta function and the generalized Riemann zeta function 

*b[0] =| [2C(-3) - C(-3, \) - C(-3, ~\) - \ (2C(-1) - C(-l, \) ~ C(-l, I)) 



IQR 

(3.35) 



This can be identified as 



3 1 11 17 , 

I6 = 6 X 240 + I20 +4 X 960- (M6) 
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the well-known sum of the zero-point energies of scalars, vectors and fermions on a 
three-sphere appropriate to the M = 4 theory. 



More interesting is the behaviour for generic nonzero (p a ij wherein we need the 
analytic continuation of the Epstein series which is known to be (see for example [f22 
or @) 



£(z,8;<pR) 



n(<pR) 



l-2z 



r z) 



r(,-l/2) + 4gc OS (ng) {7tn(pR)1/2 _ z \ ■ (3.37) 



The singularities corresponding to the individually divergent sums we encountered 
earlier appear as the poles of the gamma function T(z — 1/2) for z = — | and — |. 
Notice that these divergent poles correspond to putative mass and coupling constant 
renormalizations respectively since they are coefficients of <£> 1_2z . However, the magic 
of the M = 4 theory (whose supersymmetry is recovered in the limit (3 — > oo) ensures 



that these poles cancel against each other in (|3.34|) to leave a finite result. Hence no 
subtraction or renormalization is necessary, reflecting the UV finiteness of the theory. 



The resulting expression for the zero temperature effective potential is then 

N 

FoM=Y,C(<& j ), (3.38) 

ij=l 

with 

cV) = 

12((pR) 2 ^ K 2 (2(2n - 1)tt <pR) 2cpR (4(<pR) 2 + 1) A #i(2(2n - l)n<pR) 
tt 2 R ^ (2n - l) 2 + nR ^ 2n-l ' 

n=l v ' n=l 

(3.39) 

which manifestly vanishes exponentially for large if using standard asymptotic prop- 
erties of the Bessel functions. The reason we expect these one loop effects to vanish 
at large tpR is essentially a decoupling argument. The masses of the modes being in- 
tegrated out increase with increasing (p and hence their contributions are suppressed 
at large <p. 

For small (p, the expansion of C((p) is poorly behaved; nevertheless its asymptotic 
limit ipR — > can be extracted: 

= ( 2n -l)4 + (2n - If = 16R ' (3 ' 40) 

n=l v ' n=l v ' 

coinciding with the <p = result above. 
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Figure 4: The one-loop correction at zero temperature on S 3 . 



A numerical approximation of the zero temperature effective potential is plotted 
in Figure The full potential (along the simultaneously diagonal scalar directions) 
is obtained by adding to it, the tree level conformal coupling to the curvature. 



3.3 High temperature effective potential 



We now turn to the high temperature behaviour of the effective potential. As is 
well-known, the high temperature behaviour of gauge theories is subtle even at weak 
gauge coupling requiring a detailed understanding of the relevant length scales in 
question. 

In the high temperature limit for the M = 4 theory, the temperature is much 
larger than the inverse radius of the S" 3 

T » i . (3.41) 

Naively this is like a flat space limit. However the situation is more subtle since the 
theory has an infrared cutoff corresponding to the finite size of the three-sphere and 
what we really need to know is how the combination TR scales with the weak "t 
Hooft coupling A = g\ M N . We will come back to this issue later. For now, we note 
that a weakly coupled SU(N) Yang-Mills plasma has a hierarchy of length scales (see 
e.g fj4], Pql ), chief among these being the temperature T, followed by the perturbative 
electric screening scale or Debye length (sfXT)^ 1 and the non-perturbative magnetic 
screening length (AT) -1 . Each length scale is described by an appropriate effective 
theory obtained by integrating out higher momentum modes. 

Let us first take the naive high temperature limit (TR ^> 1) of our one loop 
effective potential ( |3.25| ). In this limit, all the mode sums over i can be approxi- 
mated by momentum integrals over a continuous variable x = nlj (TR). In addition, 
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the zero temperature contributions to the one loop potential are sub-leading and the 
high temperature effective action takes the form, 



Seff[<M] 

2(3R 3 tt 2 
where 



N 



Jr JU 1=1 n=l ij=l 



cos(n6'j 



6 \ j 

v a=l 



(3.42) 
(3.43) 



IVoiil = 

\a=l 

We have expressed our result in terms of the function 

H (y) = ~l dxx 2 e-^^ = -^{yK (y) + 2K 1 (y) + yK 2 (y)) , (3.44) 

which is plotted in ([5|). Notice that it is a negative definite, even, monotonically 

H(y) 




Figure 5: The function H(y). 
increasing function of \y\. For small y the function TC(y) approaches a constant as 



]imH(y) = -2 + ^- + 0(y 3 ) 
y~>o 2 

while as y — > oo, TC{y) approaches zero exponentially as 



hm H{y) 



-J-y^e-y 



(3.45) 



(3.46) 



The limit we have taken yields an i?-independent expression for the one loop 
term in (|3.42|) . Consequently we obtain the flat space result for a combined one 
loop potential for the scalars and the Polyakov loop in thermal Af = 4 theory. 6 

6 Similar expressions were obtained by Jl0| but for the case where the Polyakov loop variables 
are set to zero, corresponding to the deconfined phase, along with the inclusion of a finite chemical 
potential. 
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The only trace of the curved background is in the tree level mass term originating 
from coupling to the background curvature. Strictly speaking we should include an 
infrared cutoff in the momentum integrals reflecting the compactness of the three- 
sphere background, although it will not alter our result in the regime of interest. The 
regime of validity of the above potential will be discussed in more detail below. 

We will now explore some general properties of the high temperature potential 
obtained above. The one loop correction is only a function of the differences ip ai j = 
Vai — faji an d the eigenvalues <p a i each experience a quadratic classical potential due 
to the mass. Therefore for U(N) gauge group we can consistently set to zero the 
expectation value of the scalar belonging to the diagonal U(l) multiplet of the U(N) 

N 

J>« = 0. (3.47) 

i=l 

Alternatively, this is automatically true if we work with SU(N) gauge group from 
the start. Using this, and defining the dimensionless field 

Vaij = ^ — -, (3.48) 

the effective action assumes the form 

S'eff = 



N 

2(TR) 3 n 2 Y, 

ij=l 



R~ 2 1 i 

Xfittii + ~2 E( 1 -(- 1 ) n )^ n I^D COS ( n %) 
n=l 

-.2 



(3.49) 



where the 't Hooft coupling A = 9ym-N makes a natural appearance. As before we 
have used the shorthand \<f a ij\ f° r (52a ( Paij)* ■ Note the emergence of the dimension- 
less ratio involving the Debye mass scale 

R 

(3.50) 



(v/AT)- 1 ' 

which governs the size of the quantum correction relative to the classical term. The 
classical piece will dominate in the high temperature limit provided 

1 < TR < , (3.51) 
V A 

and perturbation theory in A will be valid in this regime. In other words, for pertur- 
bation theory in A to make sense, the size of the three-sphere must not exceed the 
Debye screening length scale. But since we are at arbitrarily weak coupling, we can 
always choose high enough temperature while satisfying this requirement so that the 
discrete mode sums can be effectively replaced with integrals. 
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Note also that since one loop contributions are planar, the quantum corrections 
we have computed automatically scale as N 2 in the 't Hooft large N limit, due to 
the double summation over i and j indices. 

3.3.1 A global minimum and big AdS black hole 

The function Ti(y) is negative for all values of y. Hence the absolute minimum 
of S e fi is obtained when all ip ai vanish 7 and all the 0$ are equal. Using ( |3.22| ), 
Yli ®i = m od 2tc, we find a non-zero expectation value for the Polyakov loop: 

1 1 N 

(U) = -Tre^ a = ^$>*' = e¥fc ' k = l,2,...N. (3.52) 

3=1 

Actually, in finite volume, one must define this expectation value carefully as dis- 
cussed in || in order to avoid a vanishing result due to averaging over the Ztv vacua. 
To this end one first introduces a bias which breaks the Z/v symmetry, then takes a 
large N limit and subsequently removes the bias to obtain a non-zero value for (U). 
Alternatively, one may simply compute (|?7| 2 ). A non-zero Polyakov loop implies 
that the theory is in a high temperature deconfined phase without any VEVs for 
scalars. Indeed, this saddle point of the effective potential coincides with the high 
temperature deconfined phase of the free AT = 4 theory identified in ||. 

Since we have argued that this is a global minimum of the action, we know that 
the saddle point is stable, at least within the range of temperatures in which the 
one loop approximation is valid. Importantly for us, this deconfined phase is the 
continuation to weak coupling of the large AdS Schwarzschild black hole which is the 
gravity dual of the strongly coupled deconfined plasma phase. At high temperatures, 
this minimum will dominate the canonical ensemble. 

Let us finally check that this deconfined vacuum is indeed stable to fluctuations 
in the scalar directions. Using the behaviour of Ti(y) near the origin, we find 

TV 
ij=l 

which is clearly stable for all temperatures in the range above and formally for all T. 
What we have found here is simply the finite temperature mass renormalization or 

7 Actually, there is repulsive Vandermonde type interaction which leads to an eigenvalue spread 
of 0(A) about ip — 0. This is quantum effect and subleading compared to the effects that are of 
interest to us. 



Br 



-v. 



+ 



T 2 ip 2 -- 

rat] 



as 



(3.53) 
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the thermal mass of the scalars in the deconfined phase. As anticipated earlier the 
one loop thermal mass in the deconfined phase is given by 

m th = v/AT . (3.54) 

This stability of the deconfined plasma phase at weak coupling is also in line with 
the absolute stability of the big AdS Schwarzschild black hole at strong coupling. 



3.3.2 Unstable directions charged under R-symmetry group 



So far we have restricted ourselves to a regime of temperatures ( p.51|) wherein pertur- 



bation theory in A is valid and where the classical term dominates the first quantum 
correction. In this temperature range, for arbitrary {^} the classical mass ensures 
that small ^-fluctuations are locally stable. 

Perturbation theory at high temperature 

Importantly, the one loop result actually has a much wider range of validity. It 
is in fact also valid at temperatures where it competes with the tree level term, i. e. 

Ti?~-^=, (3.55) 

while higher order perturbative contributions remain consistently parametrically 
small. The reason for this is fairly well-known in flat space: At length scales larger 
than the Debye screening length and smaller than the magnetic screening length, 
perturbation theory remains a good description with one key difference - the expan- 
sion parameter is \f\ instead of A. Of course, the situation is simpler on S 3 with 
R ~ (VAT) -1 providing an infrared cutoff. 

In the present context we may understand this as follows. In the high temper- 
ature limit, replacing discrete mode sums with integrals, the n loop contribution to 
the effective potential in naive perturbation theory will be of the form 

wWy^M ~ A^T 4 ^) (3.56) 

simply by dimensional analysis and power counting. For n > 1 these clearly appear 
to be sub-leading to the one loop term at any temperature, including when the one 
loop effect competes with the classical piece. However, a potential subtlety arises as 
ip/T —>■ 0. In flat space, at two loops and at higher orders, in this limit one encounters 
infrared divergences which necessitate a resummation of the thermal mass of the 
excitations. In our case, while we should consider a similar resummation, an infrared 
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cutoff is already provided by the size of the three-sphere. When R ~ (v^AT)^ 1 this 
IR cutoff is the Debye scale and perturbation theory continues to be valid albeit 
in the parameter v^A as opposed to A. Operationally this is automatically achieved 
since the energies of the lowest i = bosonic and fermionic modes are \f R~ 2 + <p 2 

and \J\R- 2 + ip 2 (see Table 1). Upon inclusion of this IR cutoff in our momentum 

integrals, the two loop contribution as (p/T — > scales with an extra power of \/\ 
relative to the one loop term and is suppressed at weak 't Hooft coupling. 

In summary then, within the regime of validity of perturbation theory we are 
free to consider temperatures where the tree level and one loop terms become com- 
parable. 



Unstable directions and a critical temperature 

In the temperature range R" 1 <C T < -^R^ 1 , for any configuration of 0j's, fluctua- 
tions in the ip direction are stable since the classical mass term dominates. Within 
this regime, the global minimum of the potential occurs where the absolute value 
of the Polyakov loop is unity with all ip fluctuations having a positive mass squared 
(R~ 2 + \T 2 ). 

However, the qualitative picture of the potential changes at temperatures T > 
(y/XR) . While the global minimum continues to be at ip a i — 6i — 0, the features 
of the effective potential away from the minimum can change quite dramatically. 
Expanding the action to quadratic order in <p near (p pa 0, for a generic configuration 
of eigenvalues of the Polyakov loop, we find 



S eff (0,</?)- 

N ( Rr 2 4 00 1 

2(Ti?) V £ ^— ^ + - <p 2 aij J^—[y cos((2A; - 1)%)+ ^ ^ 



Clearly, when the radius of the three-sphere becomes comparable to the Debye screen- 
ing length, it is conceivable that values of {9i} for which the cosines in the above 
formula turn negative, can lead to instabilities in the ip direction. 

Indeed, we will find that there is a critical temperature below which there are 
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no unstable ip fluctuations for a given configuration of 

T c = ^R- 1 . (3.58) 

where is an A^-dependent constant. At this temperature, an unstable direction 
develops around the configuration where the eigenvalues are uniformly distributed 
with zero expectation value for the Polyakov loop. Above the critical temperature, 
the unstable <p mode makes an appearance for a wider range of values of the Polyakov 
loop. 

Example: SU(2) 

The situation is easiest to illustrate in the SU(2) theory where the set {6^} and 
Wai}, each have only a single independent degree of freedom since 

9 X + 9 2 = 2vr, p al + p a2 = 0. (3.59) 

The SU(2) Polyakov loop is then 

U = cos9 u 0<#i<tt. (3.60) 

Rewriting the full SU(2) effective potential (|3.49|) in terms of the variables 9\ and 
ipi we find the behaviour shown in Figure ^. 




Figure 6: The effective action is plotted for T < T c on the left and T > T c on the right, as a 
function of ip and the Polyakov loop (U) for SU(2) gauge group. 

The SU(2) theory has a critical temperature 




(3.61) 
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For temperatures in the range R^ 1 <T < T c , there is a saddle point at (U) = <p = 
with one unstable direction along the [/-axis. This saddle point is actually stable 
below the Hagedorn or deconfmement temperature on S 3 || and is the confined 
phase with zero expectation value for the Polyakov loop. 



The situation changes beyond the critical temperature T c when a second unstable 
direction appears at (U) = 0. The full effective potential ( |3.49| ) now exhibits new 
extrema or saddle points with 

/lfal ^ 2 ^0 (3.62) 



T 

as in Figure |6|. It is easy to see that these saddle points must exist. At (U) = 0, we 
have seen that there is an instability at tp\ = 0. For large values of <px, the quadratic 
term in (|3.49|) dominates since TC(if) vanishes exponentially. Hence there must be 



a point for intermediate values of (tpi/T) which is stable to fluctuations along this 
direction. Importantly, in Figure |B] we see that these points are unstable to rolling 
to the (U) = 1 vacuum. Hence these new extrema are indeed saddle points. The 
instability to fluctuations in tp around <p = persists for a finite range of values of 
the Polyakov loop (U). At very high temperatures above T c , all points with tp = 
and < (U) < ^ acquire an unstable mode along the <p axis. 

SU(N) at large N 



The above picture generalizes rather directly for SU (N) gauge group. The unstable 
extremum at ip a i = is characterized by a uniform distribution of the eigenvalues of 
the Polyakov loop, 

% = ^J + c(iV), j = 1,2, ... N (3.63) 
where c is a constant chosen so that 

N 

= mod 2tt. (3.64) 

It is a stable ground state of the theory below the Hagedorn temperature on S 3 
and leads to a vanishing Polyakov loop and infinite free energy for colored sources. 
At strong coupling (and large N) and below the deconfmement transition, this phase 
is dual to the thermal AdS space saddle point of semiclassical Euclidean gravity. 

The temperatures we are considering, T > (-s/A-R) -1 at weak coupling are far 
above the deconfmement transition and hence (U) = is only an unstable extremum. 
The question we are interested in is, whether there are new unstable directions about 
this saddle point. This may be answered by examining the quadratic form in ip ( |3.57| ). 
The associated matrix is a circulant, whose eigenvalues may be determined relatively 
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easily (see Appendix A). In the large N limit, we find that as the temperature is 
increased beyond a critical value new unstable directions start to appear, with the 
first negative eigenvalue occurring at the critical temperature 

T c = -£= R~\ (3.65) 
V2X 



As the temperature is increased beyond this, we find more negative modes, each one 
appearing when the temperature hits the critical values (Appendix A), 

TP = (21 - l)-^=R-\ 1 = 1,2,.... (3.66) 
V2A 



The number of unstable directions is (N/2 — 1) when N is even and (N — l)/2 for N 
odd, and hence is formally infinite in the large N limit in the one loop approximation. 
Interestingly, this is consistent with what is observed for the small AdS black hole in 
classical gravity We should however, bear in mind that in weakly coupled field 



theory we expect perturbation theory for the constant modes or low lying harmonics 
on £ 3 to break down at high enough temperatures where the size of the three sphere 
approaches the non-perturbative magnetic screening length R ~ (AT) -1 . Therefore 
we can trust ( |3.66|) for (21 — 1) <C A=. We also point out that the modes that we 
see becoming unstable are constant on S 3 x S 1 , the higher Kaluza-Klein harmonics 
remain massive at these temperatures and within the regime of validity of the one 
loop approximation. 

As in the £77(2) example we then find new extrema where 

6 

5>-^°- (3-67) 

o=l 

The existence of these saddle points follows from the instabilities at tp = and the 
fact that at large <p the quadratic mass term dominates the potential. This can 
be verified by numerically plotting the full expression for the potential ( |3.49| ). Our 



potential is written in terms of the £0(6)-invariant radial coordinate XLv^ai i n the 
space of the six scalars. In a non-compact space, a non-zero value for this coordinate 
at the new saddle points would be interpreted as breaking the global R-symmetry of 
the Af = 4 theory to £0(5) 

£0(6) — > £0(5). (3.68) 

How can we understand this in the context of the field theory on a compact space, 
namely the three-sphere, where we do not expect to see symmetry breaking. 
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Symmetry breaking at large N 




Symmetry breaking is not normally possible in a finite volume since the wave func- 
tion spreads over the vacuum manifold or points in field space which are related by 
the symmetry. This means that we must average over the orbit of the symmetry 
group, or include in the partition function, an integral over such states. Concretely, 
introducing the notation 



1,2,... N (3.69) 



for the radial mode of scalar eigenvalues, if pi is allowed a non-vanishing value, then 
we must integrate over the corresponding orbit SO(6)/SO(5) ~ S 5 . Formally then, 
we must write the partition function as 

N 

Z= X\{d Pj p)) (cM 5 f e-^K (3.70) 

J 3=1 

Actually the above is not completely correct since the orbit of the global symmetry 
group is the symmetric product space, Sym N (SO(6)/ SO(5)) ~ Sym Ar (5' 5 ), due to 
the action of the Weyl subgroup of SU(N). However, what is important for us is 
that the space is 5iV-dimensional and the full partition function is 

r N 

Z~ Yl(dpj) e - 5 ffH+ 5 £f =1 in ft. (3.71) 

J 3=1 

However S e s scales as N 2 in the large N limit while the measure from averaging over 
the orbit of the symmetry group is of 0(N 1 ). Hence, in the large N limit we may 
neglect the contribution from the latter to compute the saddle points of the effective 
action. Note that taking the volume of the space to infinity would also have a similar 
effect since S e s scales extensively with the volume of the three-sphere. 

Symmetry breaking may now be defined in the large N theory in the conventional 
way. At finite N we introduce a source which breaks the SO (6) symmetry, 

5S = N\i ip ai . (3.72) 

Taking iV to infinity first, followed by e to zero permits breaking of the 50(6) 
symmetry. 8 



8 Global (chiral) symmetry breaking in large N gauge theory on a finite volume space has been 
studied in a different context in Bfl] 



31 



As in the SU(2) example, the unstable directions persist as we move away from 
(U) = 0, the uniform distribution for the 0j. Indeed, as some of the 9i start to clump 
together, the cosines in fl3.57| ) cease to be negative and some of the negative modes 
start disappearing. This is true in particular when the distribution develops a gap, 
an extreme example of which is the delta function distribution associated to the big 
AdS black hole which is the global minimum of the action. 



3.3.3 The gravity interpretation 



We now compare and contrast the weak coupling, high temperature picture obtained 
above of large N, M — 4 SYM on S 3 x S 1 with the known strongly coupled dual, 
namely gravity /string theory on AdSs x S 5 . 

The free M = 4 theory on S 3 analyzed in || exhibits two different phases at 
large N, a confined phase and a deconfined phase separated by a first order phase 
transition at T = T# = l/(i?ln(7 — 4y / 3)~ 1 ). It is believed that these two phases 
at A = are the infinitely curved, stringy versions of thermal AdS space and the 
big AdS Schwarzschild black hole which are separated by a first order Hawking-Page 
transition. 



In the free (A = 0) Yang-Mills theory, above Tjj the thermal AdS extremum is 
unstable. Our field theory results above show that the weakly interacting large N 
field theory develops additional (unstable) saddle points above a critical temperature 

iT 1 . (3.73) 



2A 

These configurations with scalar expectation values only preserve an SO (5) subgroup 
of the SO (6) R-symmetry of Af = 4 theory. We stress that these saddle points are 
not themselves stable, but have lower actions than the "thermal AdS" configuration 
(characterized by (U) = ipi a = 0) which develops unstable directions towards these 
new saddle points at high temperature. 



We interpret these new unstable directions as the infinite curvature or string scale 
manifestations of the Gregory-Laflamme instability of the small AdS Schwarzschild 
black hole, a phenomenon which is observed at weak curvature or strong 't Hooft 
coupling (A — > oo). As described in section 2, in semiclassical gravity this instability 
is a dynamical instability at the small AdS black hole saddle point which is already 
thermodynamically unstable. The dynamical instability triggers the localization of 
the small black hole on the S 5 breaking the SO (6) isometry to SO (5). 
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There is of course a crucial distinction between the weak and strong coupling 
regimes. The dynamical instability we see occurs at temperatures T > 0(R~ X / VA) 
which is far above the Hagedorn temperature. At these temperatures the saddle 
point corresponding to the small black hole has already disappeared |J or merged 
with the thermal AdS configuration ||. Hence it is the (U) = extremum, associated 
to thermal AdS at strong coupling, which becomes "unstable to localization on S 5 " . 

It is easy to see how the the strong and weak coupling pictures might match 
on smoothly. As the 't Hooft coupling is increased, the critical temperature ( |3.65| ) 
decreases at weak coupling. It is conceivable that as we move towards strong coupling 
the critical temperature T c at which the new instabilities are triggered, eventually 
becomes comparable to and lower than the Hagedorn temperature where the small 
black hole exists as a saddle point of the action. Hence for temperatures in the range 
T c < T < Th, the small black hole saddle point can be expected to be unstable to 
new saddle points which break the SO(6)r down to SO(5)r. As we have already 
mentioned in the introduction, there is also the possibility that at relatively small 't 
Hooft couplings with dual gravity curvatures approaching the string scale, the small 
black hole likely makes a transition to a highly excited state of strings before any 
Gregory-Laflamme like instabilities can kick in. 



4. Summary and conclusions 



In this paper we have performed an explicit computation of a joint one loop effective 
action for the eigenvalues of constant scalar fields and the Polyakov loop in SU (N) 
M = 4 SYM on S 3 x S 1 . Our perturbative computation is valid for temperatures 
< T <C jR^ 1 and in particular for temperatures where the size of the three sphere 
approaches the Debye screening length of the Yang-Mills plasma R ~ (VAT) -1 . 

We find at large N, above a critical temperature 

T c =^= R- 1 (4.1) 
V2X 



the effective potential exhibits new extrema with non-zero values for the M = 4 
scalar fields. The emergence of these saddle points is accompanied by new unstable 
directions at the thermal AdS extremum which was characterized by a uniform distri- 
bution of Polyakov loop eigenvalues. The latter is already unstable at temperatures 
below T c but new instabilities along the scalar directions appear at T c . We believe 
these unstable directions to be the small 't Hooft coupling manifestations of a well 
known dynamical instability in the gravity dual of the strongly coupled gauge theory 
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- the Gregory-Laflamme localization instability of the small black hole in AdSs x S 5 . 
Note that the weak coupling instability we have found is not associated to a small 
black hole saddle point since the latter does not exist at T c or has merged with the 
thermal AdS saddle point at the Hagedorn temperature far below T c . 

Many open questions remain, chief among these being the physical meaning of 
the instabilities of the thermal AdS saddle point at these high temperatures, and 
how they turn into dynamical instabilities of the small black hole at intermediate or 
strong coupling. Possible answers to this question have been presented in the paper 
and are summarized in a phase plot in Figure [l|. 

We also find that whenever the temperature reaches a critical value TP = 
(21 — 1)T C for I = 1,2,3,... a new negative mode shows up. A similar infinite 
set of negative modes is also expected from the classical gravity analysis of [[RJ. It 



would be interesting to understand the similarities between the two pictures better, 
particularly whether the linear dependence of TP on I can be understood at all from 
the gravity perspective, even though the latter is only valid at strong 't Hooft cou- 
pling. It is somewhat tantalizing that the values of Tc for the first five negative 



modes obtained numerically in [If:] appear to exhibit this linear behaviour. 



A natural extension of the results of this paper is to include a finite chemical po- 
tential for an 5*0(2) subgroup of the SO (6) R-symmetry. The corresponding gravity 
dual involves charged black holes. These also exhibit Gregory-Laflamme instabilities 
and we explore them at weak coupling in the holographic dual in |27J . The weak cou- 
pling analysis appears to exhibit stable saddle points with global symmetry breaking 
at finite chemical potential. The phase structure for R-charged black holes was dis- 
cussed in P?| |3D| and the dual field theory phase diagram at weak coupling has 
been explored recently in |T0| . 



The effective potential we have computed is a static effective potential since we 
calculate it in the Euclidean theory on S 3 x S 1 . It is difficult to interpret it as an ef- 
fective potential in the Lorentzian field theory on S 3 x K. since it is not clear what the 
Polyakov loops mean in the Lorentzian theory. However, there is a different analytic 
continuation to Lorentzian signatures where our effective potential may be usefully 
interpreted. This involves analytically continuing one of the angular directions on 
the three-sphere to Lorentzian signature which maps S 3 x S 1 to dS3 x S l , leaving the 
thermal circle as a spatial circle. As discussed originally in [|3l| this is the boundary 
of a "topological black hole" constructed as a certain orbifold of AdSs |H| . The 
topological black hole decays to an AdS bubble of nothing [ 33 1 where the bounce 



9 The authors of jl6| give the values of r + for each of these modes, which can be easily converted 
to a temperature. 
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solution is the Euclidean small AdS black hole which is unstable to localization on 
the S 5 . As pointed out by the authors of pi], this should be described in a dual 



field theory via an effective potential that allows rolling in a direction breaking the 
5*0(6) symmetry. The methods presented in our paper and our effective potential 
naturally allow such a holographic description of the process of gravitational decay 
via the bubble of nothing, at least in an adiabatic approximation. 
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Appendix A: Unstable directions of the effective potential 



In this appendix, we perform the small fluctuation analysis of the effective po- 
tential expanded about ip ~ for a uniformly distributed configuration of Polyakov 
loop eigenvalues ( |3.57| ). This will tell us about instabilities around the thermal AdS 



saddle point function of temperature. 

We are interested in the eigenvalues of the quadratic form in <fi<fj appearing in 
Q3.57 ) which can be rewritten as (neglecting overall multiplicative constants) 



N 



ij=l 



R~ 



1 1 



(N-l)-- + -5 NMd 
R~ 2 4 ^ 1 

xt 2 + 7T 2 (2k - 1; 



2tt 



2 cos y — (2fc-l)(i- j) 



(A.l) 



where we have used dij = 2ir(i — j)/N at the thermal AdS saddle point. 

The quadratic form above is an example of an N x N circulant matrix i.e. a 
matrix each of whose rows (or columns) can be obtained by cyclic permutations of 
the elements of one particular row (or column). The eigenvalues of a circulant 
matrix with a row of the form (x%, X2, ■ ■ ■ x^) are given by 

N 

\ e = J2x m e 2 -*r {m - 1)e £ = 0,1,2,... N-l. (A.2) 

m=l 
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Applying this formula to the quadratic form ( |A.1|) , we can find the eigenvalues 
easily. The expressions simplify at large N and we quote the results below 



A n = 



\2£-l 



A r 



Br 



XT 2 



71' 



(A.3) 



with 



1,2,.. 



It is clear that each of the eigenvalues with odd labels becomes negative at a 
critical temperature 

(A.4) 



2A 



and the first instability occurs at 



TT 



2A 



R 



-1 



(A.5) 
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